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Abstract of a Master's Thesis at the University of Miami. Thesis 
supervised by Professor Christopher N.K. Mcoers. 

Nine consecutive days of free-fall, STD data are ana'^yzed in a 
study of the baroclinic structure of the Florida Current. Mean values 
averaged over the nine days are used to reduce tidal aliasing. 

A southward flow is confirmed during the nine day period. Based 
on this data and previous works, the southward flow appears to be of 
a transient nature. 

Com.parison of the directly measured and geos trophica lly computed 
current and transport indicates that the Florida Current is essentially 
in geos trophic balance. 
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I. INTRODUCTION 


There have been several previous attempts to determine the 
validity of the geostrophic approximation in the Florida Current*, 

This thesis is a further examination of the baroclinic structure 
of the Florida Current using free-fall, STD instrument data. For 
a period of nine consecutive days in May-June 1969 direct measurements 
of the transport versus depth of the Florida Current between Miami 
and Bimini were made by Dr*> William S. Richardson (NOVA University) 
using free-fall instruments* At the majority of the stations occupied, 
one probe was equipped with the self-contained STD instrument* 

A comparison is made of the average velocity structure determined 
by differentiating the mean transport versus depth curves and of the 
structure determined by geostrophic calculations based on the mean 
density field. To first order, geostrophic equilibrium, i.e* a balance 
between the Coriolis and pressure gradient forces, holds for most 
large-scale oceanic flows* 

The first comparison of the observed and computed velocity fields 
in the Florida Current was made by Wust (1924) in which he used the 
direct measurements of Pillsbury (1890)* Considering that the geostrophic 
calculations were based on a density field determined from three 
independent sources, the agreement between the measured and computed 
velocity fields was surprisingly good* The direct measurements of the 
current field by free-fall instruments (Richardson and Schmitz, 1965) 
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have led to more recent comparisons. Broida (1966) used a quasi¬ 
synoptic density field determined from hydrographic stations, and his 
computations showed a biaxial structure in the Florida Current while 
the direct m.easurements indicated a single axis. The discrepancy was 
attributed to aliasing of the hydrographic data by internal tides. 

A time averaged comparison of data taken during the summer months of 
1965 and 1966 was made by O'Brien (1967). Using a mean T-S correlation, 
the density field was determined from the observed temperature only. 
0*Brien*s analysis confirmed the validity of the geostrophic approx¬ 
imation in the Florida Current. 

The free-fall, STD instrument provides for the first time the 
simultaneous measurement of the transport of all three of the density 
parameters: salinity, temperature, and pressure. By time averaging 
this synoptic data, tidal effects are further reduced. Thus, a rare 
opportunity exists to compare accurately the directly measured and 
the geostrophically computed velocity fields. 

With an input of geographical station locations, water depths, 
and the values of observed velocity and density versus depth, the 
output of the CHARSECT computer program (Mooers, 1970) includes the 
following baroclinic parameters: 

(1) the thermal wind ratio, which described the degree to which 
a flow is geostrophic, 

(2) the baroclinic stability parameter, which tests the criti¬ 
cality of the isopycnal slopes, and 

(3) the Richardson number, which describes the dynamic stability 


of the flow. 
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IIo METHODS 

A. DATA ACQUISITION 

The free-fall technique yields volume transport per unit width 
versus depth and surface velocity data. When the free-fall instrument 
is equipped with the STD package, a continuous trace of salinity and 
temperature versus depth is also available. The free-fall technique 
employs weighted instruments that fall (attaining their terminal velocity 
of 2 m/sec within a few meters after release) to a pre-selected depth 
where the ballast weights are released, and the instrument returns to 
the surface under its ow^n buoyancy (attaining their terminal velocity 
again within a few meters of ballast release). The precise recording of 
time and position of release and recovery provide the information 
necessary for measuring the depth-dependent transport and the surface 
velocity. At each station one drop is made to the bottom, and one to 
three drops are made to pre-selected depths. Since it is necessary to 
determine the horizontal deflection of the free-fall instrument, the 
navigational system is the controlling factor. The system used is 
Hifix (Decca Navigational System), where the master station is located 
on the vessel and the two slave stations on the western side of the 
Florida Straits. The range of the system is approximately 250 km with 
a precision of +1 meter on the western side of the Straits and +2 meters 
on the eastern side. A sire 11, high speed vessel is used in conducting 
the measurements. 

The data for this thesis were obtained from observations made over 
a nine-day period from 27 May 1969 through 4 June 1969 by Dr, W. S. 
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Richardson (NOVA University). A series of thirteen stations vrere 
occupied on a section from Miami to Bimini (Figure 1). Despite the 
short crossing time (7 to 8 hours) for a fransect as compared to a 
typical hydrographic transect (20 to 24 hours), the distorting 
influence of tidal motions remains an important factor. In an effort 
to reduce this tidal aliasing, the station times were varied as 
practicable by starting transects at Miami and Bimini on alternate 
days. Tables I and II describe the free-fall instrument data that 
were used. Those observations that included an STD drop are indicated 
in Table II- 

B. DATA ANALYSIS 

The first step in determining the density field was the digitizing 
at ten meter intervals of salinity and temperature values from each 
free-fall STD trace. At each station the mean values of salinity and 
temperature at these intervals were computed for the nine day period. 
These mean values were used as the input to the standard hydrographic 
computer program.^ The output of this program included sigma-t, 
specific volume anomaly, and dynamic depth. The cross stream sections 
of mean temperature, salinity, and (Figures 2, 3 and 4) are plotted 
with the observed current axis superimposed. 

Mean values of observed transport versus depth, rather than 
daily values, were used to compute observed velocity.^ With the volume 

^Mean value is defined as the nine day average. 

^An observed quantity is either a directly measured quantity (such 
as surface velocity) or a quantity (such as sub-surface velocity) 
directly computed from the free-fall measurement of volume transport. 
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Figure 1. Station locations. 
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TABLE I 


Station Number 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

Distance from 
reference point 
(km) 

10 

15 

20 

25 

30 

35 

45 

55 

65 

70 

75 

80 

83 

Maximum depth 
(m) 

135 

285 

345 

350 

340 

665 

765 

795 

830 

790 

735 

580 

395 

Average number 
of drops per 
station occupation 

1 

3 

4 

4 

4 

4 

4 

4 

4 

4 

4 

3 

3 

Total number of 
station occupations 

8 

9 

9 

8 

9 

9 

9 

9 

8 

8 

9 

6 

5 

Total number of 
observations with 

STD 

5 

6 

6 

6 

7 

7 

7 

6 

5 

5 

4 

3 

4 


Free-fall Instrument Data : Miami-Bimini 
Dates : 27 May 1969 through 4 June 1969 


Reference Point ; Lat. 25°44.5'N 
Long. 80008.8'W 
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STD 

STD 

Std 

STD 

STD 

stO 

STD 


STD 

STD ■ 



Time oi occupation 

0745 

0812 

0900 

0930 

1020 

1115 

1215 

1255 


1^40 

1530 1 


27 May 

Depth of probes 

122 


103 

104 ' 

102 

173 

173 

171 


340 

170 1 






176 

163 

175 1 

, 224 

361 

360 

387 


608 

332 ! 






283 

288 

323 1 

269 

448 

588 ' 

795 


780 

532 1 
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330 
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645 

755 
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STD 
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STD 

STD 

STD 

STD 

STD 

STu 


STD 

STD 


Time oi* occupation 


1428 

1355 


1320 

1240 

1120 i 

1045 

0940 

0903 

0810 1 

0750 

0655 

.28 May 

Depth of probes 


97 

86 


135 

239 1 

1 

180 

177 

169 

1 

110 

! 

176 1 

84 

86 




154 

94 


166 

3-r9 1 

341 

346 

338 

357 

354 

320 

181 




270 

329 


335 

437 

423 

563 

4 l4 

644 

512 

578 

395 

i 
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721 
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STD 
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1 

Time of occupation 

0806 

0850 

0900 

0942 

1018 

1048 1 

1150 

1245 

1340 

1442 

1512 

1600 

1630 

I29 May 

Depth of probes 

120 

88 

82 

88 

91 

173 

176 

173 

172 

165 

191 

99 

94 

1 



183 

167 

173 

164 

393 

364 

373 

329 

385 

380 1 

315 

172 




280 
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328 

463 

605 

539 

515 

563 

705 

570 

388 
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STD 
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STD 
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STD 

STD 

STD 

STD 

STD 


Time of occupation 

1348 

1324 

1248 

1218 

1148 

1112 

1018 

0918 

0830 

0748 

0712 

0636 

0612 

130 May 

Depth of probes 

118 

88 

92 

95 

95 

181 

177 

231 

233 

177 

170 

96 

98 




183 

171 

225 

175 

372 

351 

357 

344 

354 

385 

327 

177 
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265 
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494 

598 
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525 

523 
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96 

95 

93 

92 
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90 

91 
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235 

185 

134 
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208 
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372 

133 
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316 
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474 

573 
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84 

93 

145 

96 

93 

96 

81 

240 
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171 

178 

243 

167 

161 

187 

144 

169 

288 
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428 

510 

488 
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182 

367 

580 
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85 
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Figure 2. Mean temperature (C^) section 
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Figure 3. Mean salinity (%q) section 
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Figure 4 


Mean sigma-t section 
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transport available at several depths, a transport versus depth curve 
was drawn for each occupation of each station. Values were read from 
each curve at 50 meter intervals; these values were then averaged to 
obtain the mean volume transport curves. By differentiation of these 
curves, the mean velocity versus depth profiles were obtained. 

Using the method of least squares, second, third and fourth order 
polynomials w’ere fitted to these transport curves. Two constraints 
were placed on the polynomials: 

1. The transport was forced to equal zero at zero depth (sea 
surface, 

2. At zero depth, the derivative of the polynomials was forced 
to equal the observed mean surface velocity. 

An examination of the various orders of polynomials showed that the 
third order generally had the best fit and resulted in the most realistic 
profile. The third order polynomial was also used as the best fit in a 
similar treatment of free-fall data (Richardson, Schmitz and Niiler, 
1969). Comparing the velocity versus depth profiles at stations 5 and 
7 for second, third and fourth order polynomials (Figure 5), it appears 
that the third order polynomial yields the most realistic profile. 

The methods described above for determining a mean transport curve, 
and, subsequently, the velocity versus depth profiles, w’ere not the 
only ones attempted. After several tries, it was found that 50 meters 
was the minimum increment that could be used for fitting a polynomial 
to a mean transport curve. 

Another method tried was a least squares polynomial fit (with the 
same tw’o constraints) to the raw transport data, i.e. the actual values 
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Figure 5. Observed velocity profile for stations 5 and 
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of transport available at several depths,, This method resulted in 
less realistic velocity profiles, particularly in stations such as 
11, 12 and 13 where fewer observations w'ere made. 

At a sample station, direct differentiation of the mean transport 
curve resulted in a velocity profile that was comparable to the one 
obtained by differentiation of the polynomial- The method used in 
this thesis cannot be said to be the "best" method, but it does yield 
velocity profiles that are reasonable when compared with previous works. 

The resulting mean axial velocity structure across the Florida 
Straits is shown in Figure 6. The current axis of this nine day mean 
velocity field was superimposed on the temperature, salinity and 
sigma-t sections (Figures 2, 3 and 4). Of particular note is the 
southward flow near the bottom between stations 2 and 8. This southerly 
flow, which is discussed in Chapter IV, is in agreement with Wust (1924), 
bottom photos of ripplemarks by Neumann (1970), and the direct measure¬ 
ments by Duing and Johnson (1971) and Duing (1971). 

The geos trophic currents w-ere computed according to the classical 
dynamic method. As derived by Sandstrom and Helland-Hansen (1903), 
the relative velocity between two isobaric surfaces, under the 
assumption of geostrophic flovr, can be expressed as 

10 P2 P2 

V2-V1 = ^ f / 

Pi Pi 

where 

f is the vertical component of the Coriolis parameter, f = 2a)sin0 

and U) is the angular rotation of the earth and 0 is the latitude; 

Ax is the horizontal distance between stations A and B; 
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Figure 6. Mean observed axial velocity (cm/sec). 
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Oi is the specific volume, = ^35,0,p ^ ^ specific 

volume anomaly; 

and, Vi^2 pressure on isobarib surfaces 1 and 2. 

The dynamic depth, D, at isobaric surface p^^, referenced to an isobaric 
surface pg, is 

Pi 

D = / adp (II-2) 

PO 

By substitution of (II-2) into (II«1), the relative velocity between 
isobaric surfaces 1 and 2 becomes: 

V9-Vi = I^(D2 )a-(Di)a3 - C(D2)b-(Di)b^ (II-3) 

f ^ 

With the field of mass known, (II-3) can be used to compute the 
relative current field (Neumann and Pierson, 1966)• To determine the 
absolute velocity field, there must be at least one level, the reference 
level, where the absolute velocity is known* The establishment of a 
reference level is usually difficult and speculative. Computation of 
the absolute velocity field is possible if good quality direct measure¬ 
ments are available. 

The following direct and indirect methods were used to determine 
the reference level across the Florida Straits: 

(1) Determination from the observed data of the depth where the 
observed velocity equals zero. (direct) 

(2) Setting the observed and computed velocities equal at the 
surface. (direct) 

(3) Defant’s AD method. (indirect) 

The depth of no motion is the depth where the mean current component 
normal to the section is zero. The use of direct measurements to 
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determine the depth of no motion is best suited to a region where the 
mean current is strong. The observed velocity versus depth curves 
yields the depth of no motion, if one exists, for each station. Knowing 
the depth of no motion, the absolute velocity field can be computed. 

Alternatively, the directly measured surface current can be used 
to convert the relative velocity field into an absolute one. A ref¬ 
erence level results by setting the geostrophic surface velocity equal 
to the observed surface velocity. On the anticyclonic side, equating 
the two surface velocities did not indicate a depth of no motion in 
some cases, i.e. northward flow filled the entire water column.^ 

For strongly stratified fluids, one of the most reliable indirect 
methods is that proposed by Defant (1941). The Defant method consists 
of comparing relative dynamic depth differences at given isobaric 
surfaces between two stations and to determine at what level the dynamic 
depth difference is constant. This process defines a level with a 
constant horizontal pressure gradient, which is hypothesized to be the 
depth of no motion. 

The differences in the dynamic depths, computed relative to the 
sea surface, were compared for various combinations of station pairs. 

In the majority of the comparisons there was a definite layer where 
Ad was constant. Near the boundary on the anticyclonic side, several 
of the comparisons did not show a distinct level of constant Ad. 

Figure 7 shows typical Ad profiles for the cyclonic side, the middle, 
and the anticyclonic side of the Florida Current. For these profiles 

^Anticyclonic side in the Florida Current is where the relative 
vorticity is negative, i.e. the eastern side. The cyclonic side is 
where the relative vorticity is positive, i.e. the western side. 
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Figure 7. Mean dynamic depth difference (AD). 
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the layers over which Ad is constant are 60, 70 and 5 meters respec¬ 
tively. 

Using the three methods outlined above, a section showing the 
depth of no motion was computed (Figure 8). The extension of the 
reference level from m.ethod (2), Vgeos trophic = Vobserved at the 
surface, into the eastern boundary occurs where northward flow filled 
the entire water column. The doubtful areas of the reference level 
from method (3), the Defant method^ are indicated by dots superimposed 
on the solid line. It is clear that the differences between these 
m,ethods are not great. Typically it is a 25 meter difference west of 
station 8 and a 50 meter difference east of station 8. 

Different geostrophic calculations were made using various 
reference depths and various station pairs. In some cases the reference 
depth was below the maximum depth of one of the stations due to sloping 
boundaries. In these cases it was necessary to extend the dynamic 
calculations below the bottom. For these calculations, the dynamic 
depth profile was extrapolated to the necessary depth. 

Downstream velocity sections were computed using a depth of no 
motion as determined by Defant*s method (Figure 9) and by matching of 
observed and computed velocities at the surface (Figure 10). Both 
figures have the observed velocity field superim.posed, Compared to 
the geostrophic isotachs, the observed isotachs^ except near the surface, 
are deeper and skewed to the east. In both figures, the best agreement 
between observed and geostrophic occurs on the cyclonic side. 

Two additional indirect methods are shown in Figures 11 and 12. 

A velocity section was computed using the bottom as the depth of 
no m.otion. and compared with the observed field (Figure 11). The 
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Figure 8. Depth of no motion 
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Figure 9. Mean axial velocity (cm/sec), observed and geostrophic 
(depth of no motion by Defant*s method). 
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Figure 10, Mean axial velocity cm/sec;, observed and geostrophic 
(^geostrophic ~ ^observed surface). 
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Figure 12, Mean axial velocity (cm/sec), observed and geos trophic 
(depth of no motion = 400 meters). 
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biaxial structure is similar to the results of Broida (1966). Broida*s 
biaxial structure was attributed to internal tidal aliasing. In this 
case the structure is due to the use of a non-realistic depth of no 
motion on the cyclonic side. 

To illustrate the effect of a constant depth of no motion, a 
section was computed using 400 meters as the depth of no motion (Figure 
12). Again the biaxial structure is a result of an unrealistic reference 
depth. 

In each of the geos trophic sections computed, the following com¬ 
binations of station numbers w^ere used: 1-3, 2-4, 3-5, 4-6, 6-8, 7-9, 
9-11 and 10-12. All isotachs w^ere extrapolated to the boundary at the 
surface or interior. 


C o ERRORS 

The errors associated with the free-fall instrument data have two 
sources: the system errors and the possible interpretative errors 

arising from the assumptions that conditions are stationary over the 
time and space scales of a run. The system error consists of navi¬ 
gational and timing errors. The Hifix has an error of +1 meter on the 
western side of the Straits and +2 meters on the eastern side. The 
timing error is +1 second. The system error yields an error estimate 
for the surface velocity of +3% and for the transport of +5% (Richardson 
and Schmitz, 1965). For the short run time (1-10 minutes) and the short 
horizontal distance of a run (10-500 meters), experience has demonstrated 
that fluctuations large compared to the mean axial flow do not occur 


over these scales in the Florida Current. 
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A transect is completed in about 8 hours. As stated previously, 
mean values vice daily values were used to reduce the effect of tidal 
aliasing- In addition, the station times were varied as much as 
practicable to reduce tidal aliasing. Examination of the times of 
station occupation over the nine day period (Table II) reveals that 
there is little variation in those times for the stations in the 
middle of the Straits. The diurnal tidal motion may therefore have a 
greater aliasing than the semi-diurnal. Furthermore, the solar 
diurnal may have a greater aliasing effect than the lunar diurnal. 

Using the times of station occupation and values of the various 
tidal harmonic constants from Florida Current data (Smith, Zetler and 
Broida, 1969), the tidal aliasing error in the observed velocities is 
+37o for the axial velocity. The error is approximately 257o for the 
cross-stream velocity. Appendix A gives the procedure used to calculate 
these errors. 

The values of depth, temperature and salinity w’ere read directly 
from the STD trace Vvithin the limits : 
depth, +0.5 meters; 
temperature, +0.025^; and 
salinity, +0.0^o • 

The geostrophic velocity error is a sum of the relative errors in 

depth, horizontal separation, and density. The relative error in the 

density of +^0% is more than one order of magnitude greater than the 

relative errors in depth and horizontal separation. Hence, the 

relative geostrophic velocity error is +207o. 

A comparison of O' and p. at the deepest station, number 10, 

^ t '^in-situ ^ 
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showed that the approximation of the density field by resulted in 
negligible error. 

For an indication of the variability of the data, the standard 
deviations for temperature, salinity and sigma-t are plotted (Figure 13), 
The area of maximum variance in these parameters is in the pycnocline 
on the cyclonic side. 
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III, RESULTS 

A. COMPARISON OF DIRECT MEASUREMENTS 
WITH GEOSTROPHIC MEASUREMENTS 

Other than on the anticycIonic side, where the Defanc method is 
in doubt, there is a good agreement between the directly and indirectly 
determined depths on no motion (Figure 6) o From this figure, a ”best’' 
depth of no motion, called a hybrid depth of no motion, is chosen and 
used for the final comparisons of observed and computed fields. On 
the cyclonic side, the hybrid depth of no motion is the level determined 
by the Defant method. On the anticyclonic side, it is the bottom. 

Since the hybrid depth of no motion nearly coincided with the depth 
where the observed velocity equals zero, the latter could have been 
used as the best depth of no motion. The hybrid depth was chosen for 
the following reasons : 

(1) it does closely approximate the depth of no motion determined 
by the observed velocity field, and 

(2) its use sets a precedence for determining a depth of no 
motion in the Florida Straits when no observed velocity data 
are available. 

The computed surface velocity distribution using both the Defant 
and hybrid reference depths is compared to the observed surface velocity 
distribution in Figure 14. On the anticyclonic side the agreement 
between computed and observed surface velocities is poorer than on the 
cyclonic side. East of station 9, the surface velocities computed with 
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Figure 14. Mean surface velocity (cm/sec), observed and geos trophic (hybrid depth of no motion). 
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the hybrid depth shows less agreement with the observed surface 
velocities than those computed using the seemingly less reliable 
depth found by Defant's method. 

Using the hybrid depth of no motion, a comparison of the computed 
and observed velocity sections is made (Figure 15), As with the 
surface velocities, the best agreement is on the cyclonic side of the 
stream. Compared to the geos trophic velocity, the observed isotachs 
are deeper, except near the surface, and the observed current axis is 
skewed to the east. The agreement between the position of the observed 
current axis and the position of the geostrophic current axis improves 
as depth increases. 

The geostrophic, using the hybrid depth of no motion, and observed 
velocity profiles are compared in Figure 16 for each station where 
both computed and observed values were available. The geostrophic 
subsurface maximum at station 11 is similar to the observed subsurface 
maximum by Duing and Johnson (1971), The closest fit of the absolute 
values of the curves occurs on the cyclonic side. However, comparison 
of the vertical shear shows the best agreement on the anticyclonic side. 

Because of the uncertainties in determining the observed velocity 
fields by differentiation of the mean transport curves, an informative, 
if not more accurate, test of the validity of the geostrophic approx¬ 
imation is the comparison of transport curves. Figure 17 shows 
transport per unit width curves, observed and computed. A representative 
station from the cyclonic side, the middle and the anticyclonic side is 
used for the comparison. The previously established pattern of close 
agreement of absolute values on the cyclonic side and poorer agreement 
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Figure 15. Mean axial velocity (cm/sec), observed and geostrophic 
(hybrid depth of no motion). 
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on the anticyclonic side is clearly illustrated. 

The thermal wind relation derived from the assumption of geo- 
s trophic and hydros ta.tic equilibrium is 

= ^£211 
Po 

where 

V 2 is the derivative with respect to depth of the mean axial 

velocity, 

is the derivative with respect to horizontal distance of the 

mean density, 

pQ is the reference density, 
and, g is the gravitational acceleration. 

Using the observed velocity and density data, the ratio of left 
hand side to the right hand side describes whether or not the thermal 
wind relation is satisfied. If the ratio is unity, the thermal wind 
relation is satisfied. If the ratio is greater (less) than one, the 
observed shear is greater (less) than the geostrophic shear. 

From Figure 18, the areas of observed southward flow and the area 
east of station 9 are where this ratio differs the greatest from unity. 

B. FEATURES OF THE FLORIDA CURRENT 

To extract additional information about the baroclinic structure 
of the Florida Current, the free-fall data were analyzed by studying 
the following: 

(1) Baroclinic stability parameter 

(2) Richardson number 


(3) Mean T-S curves 
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Figure 18. 


Thermal wind ratio section. 



























(4) Net transport between selected isopycna Is 

(5) Total downstream transport 
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(6) Cross “Stream flow 

The results of these studies are discussed below. 

(1) The baroclinic stability parameter is the ratio of the 
slopes of an isopycnal, S, to its critical value, Sc (Mooers, 1971). 

When S > Sc, baroclinic instability, or hydrodynamic disequilibrium, 

can occur. The baroclinic stability is relatively low in the pycnocline 
on the cyclonic side (Figure 19). 

(2) The gradient Richardson number is expressed as 

Ri = — . 

(^z) 


where 

N^ = ^g is the Vaisala-Brunt frequency. 

Ri<l implies dynamic instability for the flow, and Ri>l implies 
stability. From Figure 20, the area where the dynamic stability is 
the lowest is near the bottom on the cyclonic side. The dynamic 
stability is the greatest near the surface on the anticyclonic side. 

(3) The mean T-S diagrams for stations 2, 6 and 10 (Figure 21) are 
in agreement with the T-S diagrams presented by Wennekens (1959) for 
the Florida Current. The cross channel distribution of water mass 
properties described by Wennekens is illustrated in the T«S curves. 
Station 2 is representative of the Continental Edge Water, station 6 
is representative of the Transition Zone Water, and station 10 is 
representative of the Yucatan Water. The Yucatan, or Caribbean, Water 
is identified by its well defined salinity maximum. The great reduction 
in the intensity of this salinity maximum is the conspicuous feature of 
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Figure 19. 


Baroclinic stability parameter section. 













Figure 20. Richardson number section (isolines have non-uniform 
spacing). 
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T-S curve, stations 2, 6 and 10. 


Figure 21. 
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the Edge Water. 

The water that flows through the Florida Straits originally comes, 
in large part, from the southern half of the North Equatorial Current 
and from a branch of the South Equatorial Current. This w’ater flow’s 
through the Caribbean, and then, w’ithout mixing with the w’aters 
endemic to the Gulf of Mexico, passes through the Florida Straits in 
very nearly its original state. Because this water has acquired a 
large admixture of Antarctic Intermediate Water at mid-depths from 
the South Atlantic, there is a salinity minimum (betw’een 600 and 800 m 
depth) in the w’ater exiting out of the Florida Straits (Stommel, 1966). 
This salinity minimum is present in stations 6 and 10. 

(4) The net axial transport across the Straits and between 
isopycnals was determined in the vicinity of the southw’ard flow’, A 
net transport southw’ard would favor the existence of a southward 
undercurrent rather than a large scale eddy. In all combinations of 
isopycnals, the net transport was alw’ays greater to the north. This 
result admits the possibility of a large scale eddy being the cause 
of the southw’ard flow. An example of the results are shown below. 


Isopycnal 

Interval 

Indus ive 

Stations 

Net Transport 

27.2 - 

27.3 

2 - 

9 

7.5 X 10^ m^/sec North 

27.4 - 

27.5 

6 - 

9 

3 X 10^ m^/sec North 


(5) The net downstream transport value at each station was 
integrated across the channel (Figure 22), The total mean volume 
transport of 33.4 x 10^ m /sec compares favorably with the value 
obtained by Richardson and Schmitz (1968) of 32.2 x 10^ m^/sec for data 
averaged over the period of May-June 1965. 
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Figure 22. Mean observed downstream transport (m^/sec) and surface velocity (cm/sec) with standard 
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The volume transport for each of the nine days was computed 

separately and is shown in the table below. The average of these 

b o 

nine values is 33.4 x 10^ m'^/sec. 


Da te 

Volume transport 
(x 10^ m^/sec) 

27 May 

30.3 

28 May 

33.3 

29 May 

34.7 

30 May 

33.2 

31 May 

36.8 

1 June 

31.2 

2 June 

36.5 

3 June 

31.3 

4 June 

32.9 


Continuous electrode potential measurements of transport based on 
two electrodes located at stations 2 and 3 showed good agreement with 
the free-fall measurements during the nine day period. The directly 
measured and electrode measured transport differed by less than 107o, 
but tidal aliasing of the free-fall measurements precludes any firm 
conclusion (DeFerrari, 1970). 

(6) The cross-stream data is not of sufficient quality for a 
detailed analysis of the mean cross-stream velocity structure. A 
general pattern can be ascertained by the differentiation of the mean 
transport curves. On the cyclonic side, the velocity is westward except 
for a mid-depth layer (50-70 meters thick) of eastward flow. On the 
anticycIonic side, the velocity is eastward in the upper 400-500 meters 
and westward below this depth. 









The mean cross-stream transport is shown in Figure 23, where 
positive values indicate eastward transport and negative values, 
westward transport. There is westward transport in the westward side 
of the Straits, eastward flow in the eastern side of the Straits, and, 
hence, an area of divergence near the center. The mean patterns are 
statistically significant at the 95% level. 




STATION No. I 2 ^ 4 5 6 7 8 9 10 il i2 13 

DISTA NCE ( k.t,) 1,0 . 2.0 , 3.0 , 4.0 , 5,0 , 6,0 , 7,0 , 9,0 . 
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Figure 23. Mean observed cross-stream transport (m^/sec) with standard errors 
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IV, SUMMARY AND CONCLUSIONS 

A, SUMMARY 

Nine consecutive days of free-fall, STD data, taken by Richardson 
during late May and early June of 1969, were analyzed in a study of the 
baroclinic structure of the Florida Current, The directly measured 
transport was used to obtain the mean velocity field, and the directly 
measured temperature and salinity profiles were used to obtain the mean 
density field. Mean (ioe. averaged over the nine day period) vice 
daily values were used to reduce the tidal aliasing problem. 

The mean axial velocity vras southward beneath the Florida Current 
during the nine day period. This indication of southward flow was 
consistent w’ith Hurley (1963), Neumann (1970), Duing and Johnson (1971), 
and Duing (1971). Richardson, Schmitz and Niiler (1969) did not find 
this southward flow as permanent feature when more extensive data were 
analyzed. 

For the geos trophic calculations, three methods were used to 
determine the depth of no motion: the depth where the observed 
velocity equals zero, equating the geostrophic and observed velocities 
at the surface, and Defant's method. On the cyclonic side, there was 
good agreement between the depths of no motion determined by the three 
methods, and on the anticyclonic side the agreement was not so good, 
Likevrise, on the cyclonic side Defant's method showed a definite depth 
of no motion, but on the anticyclonic side the results vrere less 
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conclusive. 

Using the hybrid depth of no motion, the geos trophic velocity 
field was computed. The comparison of the observed and computed 
velocity fields showed good agreement on the cyclonic side and much 
less agreement on the anticyclonic side. The same pattern of agreement 
between absolute values existed in the comparison of surface velocities 
and of velocity and transport profiles at several stations. The thermal 
wind equation was shown to be a good approximation, i.e. the flow was 
essentially geostrophic throughout the Straits. The greatest dis¬ 
crepancies occurred in the area of the southward flow and in the 
easternmost part of the anticyclonic zone of the Straits. The dis¬ 
crepancy may be related to the limitations of the observational and 
data analysis techniques. 

Also, the study of the standard deviation of temperature, salinity 
and showed a zone of high variability on the cyclonic side. 

The study of stability parameters showed an area of low dynamic 
and baroclinic stability on the cyclonic side. 

B. CONCLUSIONS 

(1) The validity of the free-fall, STD measurements was established 
by comparison with previous total transport values. 

(2) The Defant method was a valid method for determination of a 
depth of no motion on the cyclonic side. The bottom was a reasonable 
depth of no motion on the anticyclonic side. Thus, absolute geostrophic 
velocities can be computed independent of a direct method for depth 

of no motion determination. 

(3) The Florida Current was essentially in geostrophic balance. 
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(4) A southerly flow beneath the Florida Current was confirmed 
during the nine day period. The T-S curve of the mid-channel southward 
flow has shown the salinity minimum that is characteristic of the 
Antarctic Intermediate Water, If a steady countercurrent existed, 
different T-S curves would be likely. In the vicinity of this 
countercurrent, the net transport between isopycnals was substantially 
to the north. The lack of net southward transport, and of a distinction 
in the T-S correlations, precluded resolving whether a large eddy or 

a steady countercurrent exists on the basis of the present data. 
Richardson, Schmitz and Niiler (1968) showed that, over a longer time 
span, the north component of velocity fills the whole channel, which 
implies a transient nature for the southerly flow. The velocity 
profile analysis by Duing and Johnson (1971) and Duing (1971) gave 
definite indications of a transient southerly flow with reversals on 
a time scale as short as a day. 

Based on the present data and the previous works cited, the 
southward flow appeared to be of a transient nature and southerly 
origin. 

(5) The free-fall, STD method provided the S 3 moptic measurement 
of the velocity and density fields necessary for a study of the 
baroclinic structure of the Florida Current. Due to the tidal aliasing 
problem, time averaged mean values must be used. 

There are numerous other techniques that could have been used to 
process and analyze the free-fall, STD data. Some examples are: use 
of some other order of polynomial besides cubic; use of the spline 
technique for curve fitting rather than least squares technique; 






fitting a polynomial to the raw (i.e^ the actual values of transport 
measured) data to obtain a mean transport curve; and direct differen¬ 
tiation of the mean, or individual, transport curves rather than of 
the polynomial representative of the curves. Additionally, the same 
curve fitting technique used with the observed transport curve could 
be applied to the geos trophic transport curve. This would make the 
comparisons between the computed and observed quantities more uniform 
Because of the uncertainties in determining the **best” technique 
more emphasis in the future should be given to statistical and error 
analyses. Despite the lack of these analyses, the results of this 
thesis gave a realistic nine day mean description of the baroclinic 


structure of the Florida Current. 
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Tidal Aliasing Computations 
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The average crossing time for a transect was about 8 hours. 

Because of the distorting influence of the periodic tidal forces, it 
was necessary to work with mean values averaged over the nine days 
rather than daily values. As shown in Table II, the station times 
were varied as practicable by starting a transect at Miami and Bimini 
on alternate days in an effort to reduce the tidal aliasing. 

Because of the short record length, a tidal analysis was not 
possible. The following method was used to determine the tidal aliasing 
effect. 

For the diurnal tide, consider 



N 


and for the semidiurnal tide 



where 

cvO) is the normalized error in the mean, computed over N samples 
due to a diurnal (semidiurnal) constituent; 


6 = — 

27.3(24) 


lunar declination frequency (radians/hour); 


g _ 2 tt 

365(24) 


solar declination frequency (radians/hour) ; 


a = ilL 

25 


lunar tidal frequency (radians/hour); 



solar tidal frequency (radians/hour) ; 


N is the number of observations for a particular station; 


tj is the time of each observation for a particular station; 


and, 0 is the phase angle. 



The phase angle 0 was determined from the difference between the value 
of equilibrium argument (Vq + u) of a constituent when t = 0 and the 
epoch of a constituent (k). The values of (Vq + u) were determined 
from Schureman (1958) and values of < from Smith, Zetler and Broida 
(1969). 


Using the four constituents 0]^, M 2 and 82 ? the velocity error, 

Av, of the mean observed velocity due to the diurnal (lur.ar, solar) 
and semidiurnal (lunar, solar) tides was computed: 

^^Diurnal(lunar,solar) “ ^ ’ 

% 


and 


^^Semidiurna1(lunar,solar) ^ 

^0 


where 

are the factors described previously, 

V-p is the current speed of a tidal constituent determined from 
Smith, Zetler and Broida (1969), 

and, Vq is the ^’average" current speed of the water column estimated 
by taking % of the mean observed surface velocity. 

The values of the parameters used in the calculations and the 
resultant values of % velocity error (AV) are shown in the tables below 
Tidal constituent 0]^ M 2 S 2 

Epoch of constituent (k) 24^ 12^ 115^ 309^ 

Velocity of constituent 

(Vp) cm/sec 5.7 5.6 3.4 


2.5 






r 
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Station Number 

1 

2 

3 

4 

5 

6 

7 

Velocity of water 
column (V-p) cm/sec 

33 

68 

92 

98 

96 

89 

80 

Normalized error (cv), 
diurnal, lunar 

-.15 

-.27 

-.30 

- .27 

-.37 

-.41 

-.40 

Normalized error (cv), 
diurnal, solar 

.51 

.59 

.67 

.72 

.80 

.84 

00 

• 

Normalized error (P), 
semidiurnal, lunar 

.03 

.03 

.07 

.06 

.08 

.06 

.12 

Normalized error (3), 
semidiurnal, solar 

-.31 

-.17 

-.08 

-.06 

.17 

.32 

.40 

7o Velocity error (AV) 

4.3 

2.3 

2.4 

• 

CM 

3.3 

3.9 

4.9 

Station Number 

8 

9 

10 

11 

12 

13 


Velocity of water 
column (Vp) cm/sec 

63 

54 

51 

45 

42 

41 


Normalized error (Q^), 
diurnal, lunar 

-.36 

-.21 

-.22 

-.17 

-.17 

-.04 


Normalized error (cv), 
diurnal, solar 

.68 

.70 

.49 

.35 

.33 

.26 


Normalized error (3), 
semidiurna1, lunar 

CM 

• 

,02 

.00 

.00 

-.01 

-.03 


Normalized error (3), 
semidiurna1, solar 

.24 

.21 

-.18 

-.27 

-.32 

-.23 


7o Velocity error (Av) 

5.2 

6.6 

2.1 

0.9 

0.1 

1,9 



Hence, the error in menn observed velocity values due to tidal 
biasing is about 37o. Also, as anticipated, errors due to diurnal tides 
predominate those due to semidiurnal tides and the error is greatest 


in mid~channel. 
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and Penny Anne Stubbs. He received his elementary education in 
Statesboro Grammar School, Statesboro, Georgia, and his secondary 
education in Statesboro High School. 

In July 1958, he entered the U.S. Naval Academy, Annapolis, 

Maryland. Upon graduation in June 1962, with a B.S., he was commissioned 
an Ensign in the U.S. Navy. Subsequent naval service included two 
years in the Mine Force of the Atlantic Fleet, three years in the 
Destroyer Force of the Atlantic Fleet, and one year in the River Patrol 
Force in South Vietnam. 

He was admitted to the Graduate School of the University of Miami 
in January 1969. He was granted the degree of Master of Science in 
July 1971. 

He is married and has three children. 

Permanent address: College Boulevard 

Statesboro, Georgia 30458 
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